Evidence for BFKL and saturation dynamics from di-hadron spectra at the LHC 
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We demonstrate that rapidity separated di-hadron spectra in high multiplicity proton-proton 
collisions at the LHC can be quantitatively described by a combination of BFKL and saturation 
dynamics. Based on these results, we predict the systematics of di-hadron spectra in proton-nucleus 
collisions at the LHC. 
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Proton-proton collisions in the high energy Regge-Gribov asymptotics of QCD are described by the exchange of 
ladder like emissions of gluons. These gluons carry very small fractions x of the longitudinal momenta of the colliding 
protons but have fixed but sufficiently large transverse momenta such that the coupling as governing their emission 
is weak. In these high energy asymptotics, large logarithms of x accompany every such emission; multi-particle 
production is described by the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation P, Q which performs a resummation 
of as ln{x) terms that appear at each rung of the QCD ladder. The BFKL bremsstrahlung evolution equation in x 
can be contrasted to the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations [SHsj corresponding to QCD 
evolution with squared momentum resolution Q^. An important aspect of a quantitative description of high energy 
scattering in QCD is the regime of applicability of BFKL and DGLAP ladder resummations. 

Another important feature of QCD in Regge-Gribov asymptotics is gluon saturation [3|, which predicts that 
the non-linear dynamics of gluons slows the growth of gluon distributions at small x due to maximal phase space 
occupancies of gluons for transverse momenta below the saturation scale Qs{x) in hadron wavefunctions. This semi- 
hard scale is significantly larger at small x than the intrinsic QCD scale, thereby making feasible a weak coupling 
treatment of the non-perturbative dynamics of QCD saturation [8]. In small x asymptotics, both the growth of gluon 
distributions via BFKL evolution and the onset andproperties of the saturation regime are described by the Color 
Glass Condensate effective field theory (CGC-EFT) 

In this letter, we will demonstrate that the combination of BFKL and saturation dynamics in the leading power 
counting of the CGC-EFT provides a good description of high multiplicity di-hadron spectra in proton-proton (p-l-p) 
collisions at the LHC. These spectra were measured by the CMS collaboration jl^l in events with a high charged hadron 
multiplicity (N > 110) trigger in the rapidity range of 2 < |A?7| < 4. On the nearside (with azimuthal separations 
A(p ~ 0], the di-hadron spectra display a novel and mostly unanticipated collimation called the "ridge" . In a previous 
letter [11| , we showed that the di-hadron yield in the p+p ridge can be explained by "Glasma graphs" 12l4l4| : these 



graphs are enhanced in high multiplicity events relative to minimum bias by a^^, a factor of 10*^-10^ for typical values 
of as. We will elaborate significantly on this study here. 

On the awayside (Ai/) « vr), there is a small coUimated Glasma graph contribution, but it is dominated by the 
better known "back-to-back" QCD graphs. We will show that gluon emissions between p^'^ and p^'^ partons are 
essential and can be described by a universal BFKL Green function. Both Glasma and BFKL graphs are illustrated 
in fig. [TJ Detailed exploration of di-hadron spectra has the potential to pin down, uniquely thus far in hadron-hadron 
collisions, the underlying perturbative "QCD string" dynamics, gluon saturation, and the interplay between the two. 
In this spirit, we make predictions for di-hadron spectra in both minimum bias and high multiplicity proton-Lead 
(p-l-Pb) = 5.02 TeV collisions at the LHC. 

We first examine the Glasma Graphs for correlated two gluon production. In Ref . [iTj , we wrote down the leading 
contributions in Pt/Qsi which can be expressed in terms of the two particle momentum space rapidities i/p^q and 
transverse momenta pt, Qt as 
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where Da^^^ — *^'42(i) Qt + kT) + $^2(1) {yq, ^ ky). We extend our previous calculation by evaluating additional 
diagrams that are formally subleading in pr /Qs, as listed in appendix A of Ref. [l^. The sub-leading contributions 
discussed in that work[15j are necessary for a more quantitative description of data. The contribution of these 
additional diagrams can be expressed as 

^2 ^2'"T% ^ ^ E {Mpt,JCIt) + Mpt,J<It)) , (3) 
d^Prd^ciTdypdyg p^q^ ^\ 2 J 



3 




FIG. 1. Anatomy of di-hadron correlations. The glasma graph on the left illustrates its its schematic contribution to the double 
inclusive cross-section (dashed orange curve). On the right is the back-to-back graph and the shape of its yield (dashed blue 
curve). The grey blobs denote emissions all the way from beam rapidities to those of the triggered gluons. The solid black 
curve represents the sum of contributions from glasma and back-to-back graphs. The shaded region represents the Associated 
Yield (AY) calculated using the zero-yield-at-minimum (ZYAM) procedure. 



with[T3] Ai = 6^{pt + Qt) + 211] , such that 



^1=1 *Ai(2/p,kix)*A2(y9,PT - kij^) 

/kix 



I2 



iVp, ki_L)$A2 iVq, Pt - ki±) 



iVp, kii)$A2 iVq, Pt - kix) 



(kix ■ Pt - k^^) 
k?_L (Pt - ki_L)^ 

|ki_L X pt|^ 
k?^ (Pt - ki_L)^ ' 
(ki_L ■ Pt - k^^) |ki_L x pt| 
kf_L (Pt - ^i±f 



The other contribution, A2 in Eq. ([3]) can be expressed as 
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where kz^ = Pt ~ ~ ki_L. 

The unintegrated gluon distribution (UGD) per unit transverse area $ in these expressions is a universal quantity 
determined by solving the BaHtsky-Kovchegov (BK) equation 17, 3] as a function of the rapidity y — log(a;o/x). 
To avoid repetition, we refer the interested reader to Ref. [llj for identical details[ld| of i) the derivation of Eqs. ([T]) 
and ([3]), ii) the expression for the single inclusive distribution, iii) solutions of the BK equation, iv) a discussion of 
the non-perturbative constant C in C2 that represents the effect of soft multi-giuon interactions on the di-hadron 
spectrum. 

We now turn to the double inclusive distribution from the back-to-back QCD graphs shown in Fig. [T] At high 
energies, for a pair of hadrons having a rapidity separation Ay > l/a^, a resummation of rapidity ordered multi-gluon 
emissions is necessary, corresponding to ^-channel Pomeron exchange in the language of "Reggeon Field Theory" (20l . 
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2l| . These emissions between the two tagged partons lead to an angular decorrelation of the dihadron signal observed 



in the data. The observation of an angular decorrelation as a signal of pomeron exchange is complementary to looking 



for the growth in the dijet cross-section, as fir st p roposed by Mueller and Navelet [22[. In this framework, the double 
inclusive multiplicity can be expressed as 23], |2J| 
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where Q is the BFKL Green's function|2£ 



The <i>'s in the derivation of Eq. (O are known as impact factors and, to next-to-leading- logarithmic (NLLx) 
accuracy in x. satisfy an integral equation where the kernel is the real part of the leading-logs in x (LLx) BFKL 
kernel 23, 2^ 27 [. Since the CMS di-hadron kinematics requires very significant evolution of the $'s down from beam 
rapidities, well beyond the Mucllcr-Navclct regime, we will make the ansatz here that the <i>'s are equivalent to the 
unintegrated gluon distributions in Eqs. ([T]) and ([3]). 

As a test of this ansatz, we can "turn off' the gluon radiation between the jets by taking the as^y limit of 
the BFKL Green's function[28j and obtain for the di-hadron differential cross-section, the well known expression in 
Multi-Regge kinematics (MRK) [ii,!!!]. 
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The universal BFKL Green's function in Eq. Q resums leading logs in x, albeit with as chosen to run as the 



geometrical mean of the two momentum scales. NLLx expressions for this quantity are now known [23l |30| and 
are significantly more analytically cumbersome than Eq. (jS)). The moments Cm = (cos m(A0p_q)) ((•••) denotes the 
d'-^pjd'iq^dy dy '^sighted average over (j)p,(j3q) have been computed numerically to estimate the relative LLx and NLLx 
contributions to the double inclusive distributions for relative rapidities AY > 6 2^ 31 1. While the ratios Ci ^/Cq 
are larger at NLLx by factors of 2-3 relative to LLx, this is largely because the A(/)p ^ independent background Cq is 
significantly lower at NLLx relative to LLx. In contrast, the collimated contributions Ci and C2 at NLLx are larger 
than the LLx values by ~ 10% and 30% respectively. 

In this work, we are interested in contributions to di-hadron production that are collimated in A^p^q, above the 
^p^q pedestal. Thus while a full NLL treatment is desirable, the relatively small NLLx corrections to Ci_2 suggest 
a LLx treatment may be adequate. Further, because we include the effect of of running coupling both in the $'s 
and Q-a, piece of the NLLx contribution-the contributions neglected may be even smaller and comparable to other 
uncertainties such as choices of energy and renormalization scales. 

We now turn to a quantitative comparison of the Glasma graphs (summing contributions from Eqs. ([1]), ([3])) and 
back-to-back graphs (Eq. [5]) t o di-hadron spectra in high multiplicity collisions at the LHC. Wc first consider the 
associated yield per trigger[32| shown in the shaded region of Fig. ([T|). The detailed expressions for computing this 



yield implementing the ZYAM procedure are listed in 11|; as noted the novel feature here is the inclusion of the 



contributions in Eq. ([3]) . The results for the CMS high multiplicity (denoted here as "central p-l-p" ) events [l3| are 
shown in Fig. The agreement with data is quite good; the results in the 2 < p^'^ < 3 GeV and 3 < p?^'^ < 4 
windows have not been presented previously. While the trends are reproduced, data is slightly underestimated with 
increasing p^'® (albeit within experimental uncertainties). We note that in the absence of a Monte-Carlo generator, 
we are unable to require iVtrig > 2, as done for the data, in each transverse momentum bin. Imposing this cut would 
harden the spectrum in the correct direction. 

Inclusion of the contribution from Eq. ([3]) lowers the if-factor from K — 2.3 required previously in [11| to i^T = 1. 



One important source of uncertainty is due to the non-perturbative constant [33| (. In addition, an important caveat [34 
is that we have taken the BFKL contribution on the near-side to be independent of A(ppq for A(j)pq < 7r/2 in Figs. ([2]) 
and dS]). 

One way of addressing these theoretical uncertainties is to make predictions for the associated per trigger yield for 
p-t-Pb collisions fsH^] and compare these to forthcoming LHC data on the same. In Fig. ([3]), we show predictions for 
the nearside ridge for both minimum bias (left) and central (right) p-|-Pb collisions. All parameters are fixed to be 
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FIG. 2. The integrated associated nearside yield per trigger as a function of p!^'^ in three different p?^'^ windows for high 
multiplicity p+p coUisons. Filled symbols denote CMS data points extracted from the published pl^'^, Pt'^ di-hadron matrix p^ . 
Open symbols are our results, with dashed lines between points to guide the eye. 




FIG. 3. The associated yield per trigger as a function of p^'^ in min.bias p+Pb collisions at ^/s — 5 TeV (left) and central 
p+Pb collisions (right). Predictions are for the same pf?"^ windows as for high multiplicity p+p coUisons, and results from the 
latter are presented for comparison as dashed lines. 



the same for p+p and p+Pb. The only difference [36j is the value of the initial scale Qq in the running coupling BK 
evolution of $ for Pb nuclei. The initial saturation scale (in the fundamental representation) is taken to be three times 
that of the min. bias proton value, Qq — 0.504 GeV^, and denotes what is meant in our context by "minimum bias" 
for lead ions. A value six times larger, Qq — 1.008 GeV^, is used for "central" ions. Recall that the min. bias and 
"central" p+p values were taken to be Qq = 0.168 GeV^ and Qq = 0.672 GeV^ respectively [37,]. The key message 
from Fig. ^ is the prediction that there is a ridge in min. bias p+Pb collisions at — 5 TeV, albeit smaller than 
that in "central" p+p data. The ridge is larger in central p+Pb collisions but still smaller than central p+p data. 
As for the latter, one expects a systematic under-prediction at the higher pr,trig values due to the previously discussed 
A^trig > 2 effect. 

In Fig. (jll), we compare the "back-to-back" BFKL contribution in Eq. ([5]), which dominates the awayside A(f)pq 
distribution to data in the p^'^, p'^^^ matrixjs^. The BFKL result gives a good agreement with awayside data in the 
entire matrix. This agreement can be contrasted with the result from the MRK expression in Eq. (O, as well as the 
result of a computation in quasi-multi-Regge-kinematics (QMRK) that includes the gg — > gggg matrix element in 
this kinematics 24, 23|- In Fig. ([5]), we plot our prediction for the awayside di-hadron yield per unit A'trig ur p+Pb 
collisions at y/s — 5 TeV compared to the same in central p+p. 

It will be very interesting to learn if our predictions are consistent with observations in the forthcoming p+Pb data 
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FIG. 4. CMS data after the removal of the underlying event in each bin via the ZYAM procedure. The solid curve is the BFKL 
result added to the Glasma result with K-factors of ifsFKL = 1.1 and i^Giasma ~ 2.3. The solid (dashed) gray curves show the 
result without BFKL evolution in the MRK (QMRK) framework. 



from the LHC, and whether they can be distinguished from predictions in other related frameworks for di-hadron and 
jet spectra. With regard to the latter, it will be interesting to extend our predictions to di-jet production 4^, 41 1 



which has a larger rapidity acceptance at the LHC. Such tests of di-hadron spectra over the entire A(j)pq range have 
the potential to provide fundamental insights into the nature of QCD bremsstrahlung and gluon saturation at high 
energies. 
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FIG. 5. BFKL prediction for the away-side yield in minimum bias p+Pb collisions at = 5 TeV. The gray curves show the 
corresponding results from central p+p. The underlying event has been subtracted from each bin. 
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